Ballistic collective group delay and its Goos-Hanchen component in graphene 
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We theoretically investigate the experimental observable of the ballistic collective group delay 
(CGD) of all the particles on the Fermi surface in graphene. First, we reveal that, lateral Goos- 
Hanchen (GH) shifts along barrier interfaces contribute an inherent component in the individual 
group delay (IGD). Then, by linking the complete IGD to spin precession through a dwell time, we 
suggest that, the CGD and its GH component can be electrostatically measured by a conductance 
difference in a spin-precession experiment under weak magnetic fields. Such an approach is feasible 
for almost arbitrary Fermi energy. We also indicate that, it is a generally nonzero self- interference 
delay that relates the IGD and dwell time in graphene. 



I. INTRODUCTION 

It is well known that the process of a particle quantum 
tunneling through a barrier will take a time duration. 
Among various proposed expressions [1, 2], group delay 
[3, 4] (r g , also known as phase time in the literature) and 
dwell time [4, 5] (jd) are two well-established ones. The 
group delay is the duration between the appearing time 
of the reflection or transmission particle pulse and the 
arrival time of the incident pulse, while the dwell time is 
the time a particle lingers in the barrier region (see, Fig. 
lb). They describe the tunneling speed through a barrier 
in different aspects, and are both of paramount impor- 
tance for solid-state devices working at high frequencies 
[6] . Recently, as graphene rises as a star material in con- 
densed matter physics, extensive efforts [7-13] have been 
devoted to the investigations of group delay and/or dwell 
time in it. 

Unlike in single-electron devices, the current in bulk 
graphene devices is contributed by numerous electrons 
or holes on half of the Fermi surface. This fact implies 
that a collective group delay (CGD) of all these parti- 
cles rather than an individual group delay (IGD) of a 
single particle should be adopted to evaluate the tun- 
neling speed in graphene. The CGD can be defined as 
the summation of mode (n)-dependent IGDs weighted 
with corresponding transmission probability (T n ), i.e., 

Tg = T, n T n Tg . In this work we theoretically investi- 
gate the CGD in graphene. Our motivation is two-fold. 
(i) CGD in graphene is not directly observable. Very 
recently, an approach using a Larmor clock [5, 14] has 
been suggested to measure the transmission times r^ 
(p = 1,2, r^ 1 ) = T^ /E n T n ) at the especial Fermi en- 
ergy that aligns to the barrier charge neutrality point (or 
Dirac point) [8]. However, the tunability of the Fermi 
energy is necessary in graphene devices; while the dy- 
namics in graphene becomes totally different when the 
Fermi energy aligns to or diverges from the barrier Dirac 
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FIG. 1. (color online) (a) Sectional and (b) top schematic 
views for a particle quantum tunneling through a potential 
barrier in graphene. The red solid and blue dashed lines re- 
spectively stand for the ballistic trajectories of the A and B 
components, which are located at dzl/2/c x , a r =F l/2k x , and 
at =b l/2k x for the incident, reflected, and transmitted beams, 
respectively, (c) Experimental setup for measuring the CGD 
by spin precession. Four ferromagnetic electrodes magnetized 
respectively along the y-, y-, z(z)-, and z(z)-axis are con- 
tacted to the graphene barrier, and a weak magnetic field 
is applied in the graphene plane along the x-axis. The con- 
tact induced potential [22] can be avoided by such four points 
probe [23] and the intrinsic conductance through the barrier 
can be measured by ii4/V23|v 2 3->-o- The length of the elec- 
trodes themselves are assumed to be small compared with I 
so that the precession in them can be neglected. 



point (i.e., pseudo diffusive [8] for the former and bal- 
listic [15, 16] for the latter case). So, how to measure 
the ballistic CGD at arbitrary Fermi energy is still a ba- 
sic problem. To propose a feasible approach to measure 
the ballistic CGD at arbitrary Fermi energy is our first 
motivation, (ii) The ballistic transport in graphene is 
essentially two dimensional (2D). As a result, a lateral 
Goos-Hanchen (GH) shift [17] occurs for the reflected or 
transmitted wave packet along the corresponding inter- 



face [18-21] (see Fig. lb). As we will show, such GH 
shift makes an intrinsic contribution to the IGD. This 
fact has not been noticed in previous studies in graphene. 
It should be clarified no matter for the sake of itself or 
before we consider the measure problem of CGD. To ob- 
tain the complete expression of 2D IGD with the GH 
component contained is our second motivation. 

In this work, through the analogy between ballistic 
electrons and photons, we obtain the complete expression 
of the 2D IGD and show that, the inherent GH compo- 
nent adds an asymmetric feature to the IGD's energy 
dependence. By linking the 2D IGD to spin precession 
through the 2D dwell time, we further suggest that, for 
any Fermi energy comparable with the barrier height, the 
CGD and its GH component can be probed through con- 
ductance measurements in a spin precession experiment 
under weak magnetic fields. 



II. INHERENT GH COMPONENT IN 2D IGD 

To investigate the 2D IGD, let us consider the 2D quan- 
tum tunneling through a potential barrier in graphene 
(see, Fig. l(a)(b)). The size of the graphene sample is 
set to be smaller than the electron mean free path (£ ~ 
0.5-1/im [15, 16]) and the phase-relaxation length (L^ ~ 
3-5/xm [15]) to ensure the system stays in the ballistic 
coherent regime. A real potential occupies the region of 
< x < I and is translational invariant in the ^-direction; 
it can be induced by a top gate due to the electric field ef- 
fect [24] . The Fermi energy is E and can be controlled by 
a back gate [24] . The sample width W in the ^-direction 
is several times of I to ensure that the edge details are 
not important [25]. 

An electron with a central incident angle of a £ 
(— 7r/2,7r/2) at the Fermi surface incidents from the left 
side of the barrier. In the stationary state description, 
the spinor of electron envelope function \I/ = (ipAiipB) 
can be represented as a wave packet of a weighted su- 
perposition of plane wave spinors (each being a solu- 
tion of Dirac's equation) [18, 21]. The appearing loci 
(a?) and moments (t^ 77 ) for the A and B (77 = ±) 
component of the incident, reflection, and transmission 
(£ = i,r,i) packet peaks (see, Fig. 1(b)) are deter- 
mined by the following condition: the gradient of the 
total phase (cj)p = <j)^ + k y a7 — Erp /K) in the wave vec- 
tor (fe = k x x-\-k y y) space must vanish. Here, r = |r|e^ r 
(t = |£|e^*) and <j) r {cj) t = <j)t + k x l) are the reflection 
(transmission) coefficient and corresponding phase shift, 
respectively (see, Fig. 1(b)). This is similar to the opti- 
cal case [26] but with the extra sublattice degree of free- 
dom (77). Since k x = Ecosa/hvp and k y = Esina/hvp, 
the above condition means two independent conditions 
d(j)p jdE = and d(f)p /dk y = 0. A comparison of the 
second conditions between the reflection (transmission) 
and incident beams gives the sublattice dependent lat- 
eral GH shift [18, 21]: af = -{d<j> r /dky) E T l/k x and 
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FIG. 2. (color online) The 2D IGD, its scattering, and GH 
components (in units of the equal time to = 1/vf) as a func- 
tion of the Fermi energy at a = 20°. The magenta thin 
solid curve stands for the time associated with the Snell shift, 
rs = <?s sin a/vF ■ The parameters of the potential barrier 
(which are also used in Figs. 3-5) are l/lo = 1 and V/Eo — 3tt 
with /o and Eq = Hvf/Io being a length unit and energy unit, 
respectively. 
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- —{d(j)t/dky)E (see, Fig. 1(b)). A comparison of 
the first conditions gives the reflection and transmission 
2D IGD, 



di^ckyOj) 



dE 



(1) 



Here sublattice independent average values of GH shifts 
o-£ = —(d(j)^/dk y )E have been used to evaluate the IGD 
for the electron spinor rather than its sublattice com- 
ponents. For asymmetric barriers, there is a differ- 
ence between the reflection and transmission in both <p^ 
and <j£, so a bidirectional 2D IGD can be defined as 
T ff = Ejle| 2 r|. 

In the ID case the bidirectional IGD has an expression 
of Tg f = Y^t \i\ 2 ^d(j)^/dE^ which stems solely from the 
phase shifts introduced by scattering at the interfaces 
[27]. Comparing the two expressions, we can find that 
the differences come in two aspects between the 2D and 
ID cases. Firstly, the differential in the ID case becomes 
partial differential in the 2D one, i.e., d/dE — >• (d/dE) a . 
Secondly, besides the scattering phase shifts in the ID 
case, the GH shifts also contribute inherent phase shifts 
in the 2D case, i.e., <p^ — » <p£ + k y a^. 

To highlight the intrinsic contribution of the GH shifts, 
we rewrite the bidirectional 2D IGD as 



st 



GH 



(2) 



The first component r^ = J\ \^\ 2 h(d(j)^/dE) a comes 
from the scattering phase shifts at the interfaces, which 
is the same as the ID case; while the second component 
rf H = J2? I£| 2(J £ sina/vp results from the GH shifts, 
which is totally absent in the ID case. Note, we have 



used the relation h(dk y /dE) a = sina/vp to get the GH 
component, since the GH shifts are not explicit functions 
of the energy. We would like to stress that, although 
the GH component displays a partial differential respec- 
tive to k y at fixed E (i.e., r^ H = J\ |£| 2(7 £ sina/vp and 
<j£ = —{d(j)£ t /dky)E), this component essentially comes 
from the condition that the energy gradient of the GH 
phase shifts (i.e., d/dE) must vanish. We can understand 
the GH component as following. At time h(d(j)t/dE) on 
the wave front of the transmitted beam reaches (/, 0) and 
then propagates freely with a velocity of vp to the final 
position (l,a t ) (see, Fig. 1(b)). This step will cost a 
duration of a t sin a/vp since the wave front is perpen- 
dicular to the propagation direction. This picture also 
holds for the reflected beam. A weighted average of the 
transmission and reflection gives r^ H . 

Figure 2 shows clearly the contribution of the GH com- 
ponent to the 2D IGD. One can see that rj? 1 is sym- 
metric about the center (E/V = cos -2 a) of the trans- 
mission gap (TG) due to the symmetry of ^ about it. 
On the contrary, r^ H stemming from the quantum GH 
shifts is asymmetric about the TG's center. The quan- 
tum GH shift displays the same trend as the classical 
shift predicted by the Snell's law (as = I tan j3 with 
f3 = sin - (hvpk y /(E — V)) being the refracted angle). 
The latter is obviously asymmetric as it is negative (pos- 
itive) in the low (high) energy range (see Fig. 2). In total, 
the GH component not only quantitatively contributes a 
part of order of as sina/vp to the 2D IGD, but also qual- 
itatively results in the remarkable asymmetric feature in 
the energy dependence of the 2D IGD. 



III. MEASURING THE CGD BY SPIN 

PRECESSION-INDUCED CONDUCTANCE 

DIFFERENCE 

Having obtained the complete expression for the IGD, 
we now seek physical observable for the CGD at arbitrary 
Fermi energy. The Larmor precession of the electron spin 
in a magnetic field provides a clock for studying the elec- 
tron dynamics [5, 8, 28]. Here, we consider a configura- 
tion where the magnetic field is applied in the graphene 
plane along the x-axis (see, Fig. 1(c)). In such a con- 
figuration, the dynamical perturbation by the Lorentz 
force is avoided, and the only effect of the magnetic field 
is to cause spin precession round the x-axis. For elec- 
trons with spins initially directing the ?/-axis, a duration 
of this precession leads to different transmission proba- 
bilities in the z and z directions. We calculate the spin- 



dependent transmission probability (T z ^)y 



% 



(z)y\ 



*)in 



a pseudospin-spin direct product space (a similar method 
is used in Ref. [8]), where T z ^) y denotes the transmission 
probability for an electron incident from electrode 2 (with 
^/-directed spin) and transmitted to electrode 3 (with z 
or z-directed spin) (see, Fig. 1(c)). An explicit equality 
between the spin polarization in transmission probabil- 
ities P = (T zy — T zy )/(T zy + T zy ) and the dwell time 




FIG. 3. (color online) The dwell time and P/ujb at various B 
as a function of the Fermi energy at a = 20°. 



(rd = J \^(x)\ 2 dx/(vpcosa)) is found in the weak- field 
limit: 



r d = (P/uj b )\b- 



(3) 



where, uob = Q^bP is the Larmor frequency, g is the 
gyromagnetic factor in graphene, jjlb is the Bohr magne- 
ton, and P is the magnetic field with a reduced strength 
B = Hujb/^Eo. This equality is clearly demonstrated 
in Fig. 3. Note, this equality is restricted in symmet- 
ric structures [29]. Under such restriction Eq. (3) holds 
for arbitrary incident energies and angles and can be in- 
terpreted physically in following. Rewriting the equality 
as UBTd = P and multiplying both sides by ft/2, the 
right hand side gives the expectation value of S z of the 
transmitted electrons. Such expectation value is deter- 
mined by the product of the spin precession frequencey 
(i.e., the Larmor frequency) and the time the precession 
persists (i.e., the dwell time). This is just what the left 
hand side expresses. It is noted that Eq. (3) also holds 
in systems with parabolic dispersion relations and scale 
envelope functions [5]. 

To obtain the summation of the dwell time on half of 
the Fermi surface, we multiply Eq. (3) by (T zy + T zy )/2 
(which tends to T when B — >• 0). Then the right hand 
side of Eq. (3) can be further related to an experimen- 
tal observable, the conductance (G). This is because 
the conductance is defined as G(E) = Go^2 n T n , where 
Go = 2e 2 /h is the quantum conductance considering the 
twofold valley degeneracy. Thus, Eq. (3) can be rewrit- 
ten as 



E^ B)r - 



G zy (E) - G zy (E) | 
2ujbGq 



|s^o, 



(4) 



where G z ^ y can be probed by Im/V^Iv^-^ hi the pro- 
posed experimental setup (see, Fig. 1(c)). 

We now try to link the CGD to the above conductance 
difference by the relation between the 2D IGD and the 
2D dwell time. This relation can be obtained by making 




FIG. 4. (color online) The reduced IGD, dwell time, self- 
interference delay, and transmission probability as a function 
of the Fermi energy at a = 20°. Insert: the collective self- 
interference delay versus the Fermi energy. 



a variation of the time-independent Dirac equation at 
the barrier boundary with respect to the two individual 
variational parameters, (k x , k y ) or equivalently (E, a). 
The 2D feature of the tunneling process and the spinor 
nature of envelope functions should be taken into account 
in the variation, for which the detailed derivations can 
be found in the Appendix. We make the variation with 
respect to E and sin a, since the variation results about 
them can be expediently related to r^ (see, Eqs. (A3)- 
(A7) in the Appendix) and rf H (see, Eqs. (A8)-(A9) 
in the Appendix), respectively. The concise result for a 
rectangular barrier reads 



Td 



T q +T i: 



(5) 



where t\ = h\Re(r) cos a + Im(r) sin a] sin a/(E cos 2 a) is 
a self- interference delay (see, Eq. (A7) in the Appendix) 
stemming from the interference of the incident and re- 
flection envelope functions in front of the barrier [1, 30]. 
The correctness of this relation can be verified by numer- 
ically calculating and comparing the explicit expressions 
of rd, T g , and n. 

Fig. 4 shows the 2D IGD, dwell time, and self- 
interference delay in reduced form as a function of the 
Fermi energy at a fixed incident angle, where e = E/ Eq 
is the reduction factor. One can see that the reduced self- 
interference delay (— eri) is symmetric about the TG's 
center; it achieves the maximum at the TG's center 
and oscillates around zero outside the TG. The self- 
interference delay itself is important only in the low en- 
ergy range (diverging as E~ x when E — »> 0). It disappears 
at (anti)resonant tunneling since there is no interference 
in front of the barrier. Accordingly, the IGD almost co- 
incides with the dwell time except within the low energy 
range or around the TG's center. 

The collective self-interference delay (rf 7 ) can 

be calculated by using ^2 n T n (E)r^ n (E) = 
Mf^ 2 /2 Ti(E,a)T(E,a)d(sma), where M 




FIG. 5. (color online) The CGD, its scattering component, 
and the H-dependent conductance difference (i.e., the right 
hand of Eq. (6)) as a function of the Fermi energy. Insert: 
G zy (solid) and Gzy (dashed) could be directly measured in 
the experiment for B — 10 _1 . 



\E\W/hv F = (\E\/Eo)(W/2tt1 ) = M E M W is the 
number of the transverse modes (which should be rather 
big). As can be seen in the insert of Fig. 4, rf oscillates 
with E with the amplitude disappearing exponentially at 
a relatively high Fermi energy. Therefore, Eqs. (4) and 
(5) imply that, for arbitrary Fermi energies comparable 
with V, the CGD can be directly measured by the 
spin-precession induced conductance difference under 
weak fields 



ifm 



Gzy — Gz 

2ujbGo 



B^O- 



(6) 



Fig. 5 shows the CGD, its scattering component, and 
the magnetic field dependent conductance differences as 
a function of the Fermi energy. To clearly show the oscil- 
lation details in the low energy range, all these qualities 
are plotted for a single mode (i.e., divided by M). As 
is seen, the spin-precession induced conductance differ- 
ence can be obtained by measuring G zy and G zy sepa- 
rately (see insert in Fig. 5) and it increases for weaker 
B (see Fig. 5). At B = 10 -3 , it is already a rather 
good measurement of the CGD for E > 0.2V. Thus, we 
provide an electrostatic approach to measure the CGD 
in graphene, which is feasible for almost arbitrary Fermi 
energy. This approach looks interesting when we notice 
that the CGD describes a dynamic process while the con- 
ductance is static. 

In Fig. 5, it is also noted that, due to the GH shifts 
with different signs, the scattering component is always 
larger (smaller) than the expected value of the CGD 
when E < V (E > V). Then, a comparison between the 
experimentally observed value of the conductance differ- 
ence (the right hand side of Eq. (6)) and the theoretical 
prediction of the CGD (the left hand side of Eq. (6)) 
at Fermi energies away from the barrier Dirac point can 
be utilized to probe the inherent GH component of the 



CGD (or equivalently the intrinsic effect of the GH shifts 
on the CGD) in graphene. 



IV. CONCLUSIONS AND REMARKS 

In summary, we have obtained the complete expres- 
sion of the 2D IGD in graphene with the intrinsic contri- 
bution of the GH shifts being contained. Moreover, we 
have proposed an approach to probe the ballistic CGD 
and its inherent GH component by conductance measure- 
ments in a weak-field spin precession experiment. This 
approach is feasible for almost arbitrary Fermi energy. 
We have also indicated that, it is a generally nonzero 
self-interference delay that relates the 2D IGD and dwell 



time in graphene. The inherent GH component of the 
2D IGD should also present in other 2D ballistic coher- 
ent electronic systems, since the derivations of it do not 
depend on the specific electronic excitation of graphene. 
The feasibility of the proposed approach to measure CGD 
in these systems is an interesting question. 
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Appendix A: Derivation for the relation between 2D 
IGD and dwell time in graphene 

Let us begin with the single-particle Dirac equation 
that governs the low-energy excitation in graphene. In 
the barrier region it reads as [vpcr • p + V(x)] 1 9(x,y) = 
E*&(x,y), where the pseudospin matrix a has compo- 
nents given by Pauli's matrices and p = (p x ,p y ) is the 
momentum operator. The eigenstates ^(x,y) are two- 
component spinors with each component being the en- 
velope function at sublattice site A/B of the graphene 
sheet. Due to the translational invariance along the y- 
axis, the envelope function can be separated as \l> = 
[^a(x), ^b(%)Y z lkvV with k y = Esma/Hvp- The A and 
B components along the x-direction are related by a pair 



of coupling first-order equations 

d E-V 
-z-iI>a,b = ±k y ^ A ,B 77 ll>B,A, 

OX ihVp 



(Al) 



which imply a decoupled second-order equation for both 
the A- and 5-components 



dx 2 



2^1 



(E - Vf - E 2 sin^ a 



(hv F y 



^A,B = o. 



(A2) 



We carry out the energy- variational form and conju- 
gate form of Eq. (A2) and upon integration over the 
length of the barrier we get 



A,B 



d^JA,B 9t 
dE dx 

1 2Ecos 2 a 



■r A 



2 ^A,B 



I 

Jo 



A ' B dEdx 
2V(x) 



X — l 



x=0 



(A3) 



(hv F y 



-IV> 



A,B\ 



J dx. 



It is found that when is evaluated by the envelope func- 
tion outside (inside) the barrier, the left (right) part can 
be related to r^ (r^). However, we should note that Eq. 
(A3) is only valid inside the barrier as the spatial deriva- 
tive of ipA,B are not continuous on the potential bound- 
ary. To overcome this dilemma, we express dipA,B/dx 
inside the barrier by Eq. (Al) and their conjugate form. 
Since ipA,B themselves are continuous, envelope functions 
inside the barrier can be replaced by the corresponding 
ones outside the barrier. Then the left part of Eq. (A3) 
can be evaluated. For the A-component it reads as 



J + K + cos a(—ir + ir*)/hvp 1 
and for the £> -component, the result becomes 
J + K + cos a(ie~ i2a r - ie i2a r*)/hv F , 



(A4a) 



(A4b) 



wheTeJ = ^{[B(0)-A(0)}\r\ 2 ^ r ^[B(l)-A(l)}\t\ 2 ^} 
and K = j^{[B(0) - A(0)]|r||rf + [B(l) - ^(0]l*ll*r>- 
Here the relation of lossless barriers \t\ 2 + \r\ 2 = 1 
has been used and the notations are adopted as O' = 
(dO/dE) a , A{x) = sin a + iX(x)e za , B{x) = sin a - 
i\{x)e~ ia , and X(x) = 1 — V(x)/E, a ratio of the kinetic 
energy inside and outside the barrier. 

Since i/j^i/j = ^\^a + ^b^b, the relation for each 
spinor component should be added to get the relation 
for the spinor, which at last reads 

f Q [\(x)-sm 2 a]\^(x)\ 2 dx 
vf cos a 
=A(0)|rf fi<#. + \(l)\t\ 2 h(f)' t - i\(0)h\r\\r\ f - i\(l)h\t\\t\' 
[Re(r) cos a + Im(r) sin a] sin a dk 
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dE' 



(A5) 



where Re(r) (Im(r)) is the real (imaginary) part of r. 
We can see that, this equation is a general result that 



relates the integral of the weighted probability density 
inside the barrier (left part) and the weighted energy- 
variational behavior outside the barrier (right part). It 
is noted that, the factor X(x) has a critical role in the 
relation. 

To clearly relate the general result Eq. (A5) with both 
Td and Tg f , we consider a restricted condition that X(x) is 
a constant under the barrier (i.e., a rectangular barrier). 
Note this condition is not necessary for the common semi- 
conductors case [27], a reflection of the spinor nature of 
graphene. Under such condition, the third and fourth 
terms on the right hand side of Eq. (A5) disappear due 
to the lossless condition of the barrier |t||t| / + |r||r| / = 0, 
and Eq. (A5) can be rewritten in terms of r^ and r^, 
i.e., as a sub-relation 

Td(X — sin 2 a) = r^X + Ti cos 2 a, (A6) 

where a self-interference delay is found from the last term 
ofEq. (A5), 



h\Re(r) cos a + Im(r) sin a] sin a 
E cos 2 a 



(A7) 



Graphene is two dimensional, which means there are 
two independent variation parameters, k x and k y or E 
and sin a. We have made the variation about E, the 
variation of Eq. (A2) about sin a reads 
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Following a similar way as above, we straightforwardly 
obtain the sub-relation between r^ and r^ H under the 
same restricted condition of a constant A. The sub- 
relation reads 



Td sin a = Tq X-\- Ti(X — cos a). 



(A9) 



Making a simple addition of the two sub-relations in 



Eqs. (A6) and (A9) and taking into account r g = r ( 
rf H (see, Eq. (2) in the main body), we finally get 



st 



Td 



(A10) 



The correctness of this relation and the sub-relations can 
be verified by numerically calculating and comparing the 
explicit expressions of the five times. It may be valu- 
able to indicate that, many previous works [9-11] have 
also concerned such relation and they all have achieved 
a same conclusion that r g = Td (actually they mean 
r g t = Td smce T f H i s n °t included in their discussions). 
The results obtained here clearly indicate that, it is a self- 
interference delay that relates the group delay and dwell 
time in graphene, similar to the common semiconductor 
case [27]. 

For the normal incident or ID tunneling case (a = 0), 
T{ vanishes (see the factor sin a in Eq. (A7)), since no 



reflected portion thus no interference happens in front of (A7) also hold for the tunneling of massless Dirac parti- 
the barrier due to the Klein tunneling [31]. The rela- cles in topological surface states [32] , where the real elec- 
tion and subrelations revealed in Eqs. (A10), (A6), and tron spin rather than the sublattice structure in graphene 
(A9) and the expression for the self-interference delay Eq. provides the Dirac structure. 



